Quantum creep and quantum creep transitions in ID sine-Gordan chains 
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Discrete sine-Gordon (SG) chains are studied with path-integral molecular dynamics. Chains 
commensurate with the substrate show the transition from collective quantum creep to pinning 
at bead masses slightly larger than those predicted from the continuous SG model. Within the 
creep regime, a field-driven transition from creep to complete depinning is identified. The effects 
of disorder in the external potential on the chain's dynamics depend on the potential's roughness 
exponent H, i.e., quantum and classical fluctuations affect the current self-correlation functions 
differently for H = 1/2. 
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The discrete sine-Gordon (SG) chain, also known as 
Frenkel-Kontorova (FK) model 00, is a generic model 
for the motion of an elastic object composed of discrete 
degrees of freedom through an external potential. To 
name a few of its applications, the SG model is used 
for the description of driven charge-density waves in 
solids, coupled Josephson junctions 0, the sliding mo- 
tion of an adsobed layer of atoms over a subtrate [4( , and 
most recently electronic conductance in nanotubes 0, 
see also the reviews 0, on electronic transport in 
one-dimensional structures. A lot of attention has been 
devoted to the (discrete) classical FK model both at 
zero and finite temperatures 01 and the (continuum) 
quantum- mechanical SG model 0,0. However, less is 
known about the quantum-mechanical properties of the 
discrete quantum FK (QFK) model, in particular about 
its dynamical properties. Numerical approaches have 
lead to a clear picture, how quantum fluctuations renor- 
malize the thermal equilibrium structure, but results are 
often limited to zero (or small) external fields or to 
variational approaches While quantum Monte Carlo 
(QMC) simulations yield numerically exact results for 
static properties, they only allow one to calculate small- 
frequency dynamical properties indirectly, i.e., conclu- 
sions on the existence of a phonon gap are drawn by 
studying the temperature dependence of the internal en- 
ergy In some cases, more dynamical information 
can be withdrawn from QMC if the functional form of 
the low-energy spectrum is known 

It is well established for the SG model that the effects 
due to thermal and those due to quantum fluctuations 
differ qualitatively. Thermal fluctuations automatically 
lead to creep, while quantum fluctuations do not. At fi- 
nite T, kink anti-kink pairs will be activated and a small 
external force will eventually be able to drive the pairs 
apart, resulting in net mass transport 0. If one assumes 
random interaction between the elastic object and the 
embedding' svstcm. thermal creen is also nresent, in higher 



dimensions, as shown in a sophisticated renormalization 
group (RG) study H3- Strictly speaking this implies that 
the pinning or the static friction force F s is zero and that 
the kinetic friction force vanishes linearly with sliding 
velocity vq in the limit of v — > 0, even thoug the propor- 
tionality coefficient 7 ff = lim„ ^o -Ffc/uo may be large. 
In order for the ID, quantum sine- Gordon model to creep, 
it is not sufficient to have arbitrarily small quantum fluc- 
tuations, but the effective masses m (defined as density 
times period of substrate potential) must be less than 
a certain critical value m c (at fixed momentum cut-off 
and fixed substrate strength) [8|. For m > m c , thermal 
fluctuations and/or finite external forces are required to 
initiate mass transport. Recent RG studies 0| suggest 
that m c is also finite in higher dimensions if the elastic 
manifold is pinned through an external random potential. 

In many of the above mentioned cases, including the 
quantum SG mo del, it is necessary to go beyond standard 
RG theories |lj, Ijjj [lj| , because otherwise the effect of 
creep or alternatively the effect of a vanishing excitation 
gap at zero wave vectors might be artificially removed. 
It is thus desirable to have numerical techniques which 
allow one to (i) verify the results of RG studies, (ii) ob- 
tain results beyond continuum approximations allowing 
an accurate determination of critical values such as m c 
for discrete systems, and (iii) obtain dynamical responses 
of the system for arbitrary external forces. In this Letter, 
we first intend to establish that (adiabatic) path integral 
molecular dynamics (aPIMD) [rH Hif is a well-suited 
technique to tackle the many-body dynamics of elastic 
manifolds moving through embedding systems. We will 
then study the FK chain's dynamics as a function of the 
external field. This analysis includes new results for the 
quantum dynamics of a FK chain with various types of 
disordered substrates. 

PIMD is a method to calculate static 0] and dy- 
namic 19] properties of many-body quantum systems 
in thermal enuilibrium. Tt is based on the isomornhism 
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between the partition function of a quantum mechani- 
cal point particle and that of a classical ring polymer. 
The number of monomers per ring is called Trotter num- 
ber P. Recent progress was made concerning the dy- 
namical interpretation: It was suggested that real time 
correlation functions C (t) of observables linear in veloc- 
ity v and/or position x can be calculated exactly with 
adiabatic PIMD methods [3 HJ]. The C(i)'s follow 
from the equivalent correlation functions C c (t) defined 
for the ring polymer's centroids (center of mass), pro- 
vided that the P — 1 inner degrees of freedom are in 
full thermodynamic equilibrium for every given position 
of the centroids. Furthermore, the characteristic time 
constant of the centroids' thermostats must be large 
compared to the intrinsic relaxation time of the sys- 
tem. If these two conditions are satisfied, C(t)'s and 
C c (£)'s Fourier transforms are related through the equa- 
tion C(w) = f3huj/2(coth(f3huj/2) + l)C c (w), where /3 is 
the inverse thermal energy. In the following, only cen- 
troid spectra C c (ui) will be considered, as they reflect 
directly the density of states. 

While we refer to Ref . EH and to a future longer- version 
publication for more details of our implementation, we 
may note that we used a representation of our system 
that is diagonal in the harmonic part of the Hamilto- 
nian H. Moreover, higher-order corrections to the high- 
temperature density matrix were included in our treat- 
ment leading to P~ 4 convergence j^. We made sure 
that the relative systematic errors of the internal energy 
with respect to the classical ground state due to the use 
of finite P was less than 10" 3 . While PIMD is a finite- 
temperature method, it is possible to extract zero-T be- 
havior by analyzing the convergence of the results by de- 
creasing T, just like the thermodynamic limit can be ap- 
proximated from finite-size extrapolations. We will first 
consider the most simple QFK model, which corresponds 
to a discrete, elastic chain, which is commensurate with 
the underlying potential or substrate. The Hamiltonian 
H reads 
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where p n and x n are, respectively momentum and posi- 
tion of particle n, fc is the stiffness of the spring con- 
necting two neighbored particles, Vq is the coupling 
strength to the embedding system, and 2irb is the sub- 
strate's lattice constant. The periodic boundary condi- 
tion i n+ n ~ x n + 2TtbN makes the chain commensurate 
with the substrate. Our system of units will be defined 
by Vq, 6, h, and Boltzmann's constant fog. Unless other- 
wise noted, we vary the mass m but leave the harmonic 
intrachain coupling k = 0.1 Vo/b 2 constant. This value 
of k is much smaller than the maximun curvature of the 
potential max(d 2 V(x)) — Vo/b 2 . This challenges the con- 
tinuum approximations, which assume slow variations of 
the reduced positions x n — 2nirb with index n. 
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FIG. 1: Left: Excitation gap cu(q = 0)(m) at zero wavelength 
as a function of mass m. The arrow indicates where the con- 
tinuum sine- Gordon model shows the transition from finite to 
zero gap. Right: Phonon dispersion for two different masses 
and for the classical case. Lines are fits according to Eq. @ 
where ko is the only fit parameter (ujo = y/ ko/m). 



Classically, the dispersion relation of the chain is sim- 
ply given by 
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where q denotes the phonon's wavenumber and ko = 
ko '■= Vo/b 2 . Early calculations |23| suggested that this 
classical, zero-temperature value is renormalized to a re- 
duced effective coupling fco due to zero-point and/or ther- 
mal fluctuations. Also for the discrete, commensurate ID 
QFK model, Eq. J2J with a non-trivial value for ko turns 
out to provide an excellent approximation to our calcu- 
lated phonon spectra. The dispersion relations of our 
chain are shown in the right-hand side of Fig. ^for two 
different masses. (The dispersion relation was obtained 
by fitting Lorentzians to spectra similar to those shown 
below in Fig. [3] In the commensurate system, the lines 
are much sharper than in the disordered case.) The clas- 
sical zero-temperature dispersion relation is inserted for 
comparison. Each curve requires only one fit parameter, 
namely the value for fco ■ The dispersion relation does not 
change when particle number N and Trotter number P 
(sufficiently small but fixed (3/P) are increased. We want 
to note that all relevant imaginary-time correlation func- 
tions could be well reproduced from the centroid spectra 
C c {ui), which made us confident to use the dynamic in- 
terpretation of adiabatic PIMD. 

The phonon excitation gap apparently becomes zero 
at a value m c w 0.02 h 2 /Vob 2 , as shown in the left- 
hand side of Fig. ^ This in turn implies that sliding 
can be induced with an arbitrarily small external driv- 
ing force for m < m c . As we will show later, the sys- 
tem creeps in the zero-gap regime when subjected to 
a small external driving force. The discreteness of the 



chain alters the value of the mass m c at which the tran- 
sition from finite gap (no creep) to zero gap (creep) 
takes place. The continuum model predicts this tran- 
sition to occur at m c = 0.016 h 2 /V b 2 24]. Thus for 
k = 0.1 Vo/b 2 , the value for m c differs by about 20% 
between the discrete FK and the continuum SG model. 
This discrepancy will decrease as the spring stiffness 
within the chain increases as compared to the maximum 
curvature fco = Vo/b 2 of the embedding potential. With- 
out showing the data explicitly, we would like to com- 
ment again on the imaginary-time behavior, as this be- 
havior does not rely on the correctness of the claims made 
for adiabatic PIMD. For m > m c , we observe a finite 
plateau value in the imaginary-time correlation function 
G(An, At) = ({x„ +A n(r + At) - x n (r) - 27r&An} 2 } for 
An and At — > oo. Similar to what one would expect 
from a (quantum) Berezinski-Kosterlitz-Thouless transi- 
tion |2al26j, G(An, At) increases logarithmically both 
with An and At for m < m c . An accurate determina- 
tion of m c 's value would yet remain much more difficult 
in terms of an imaginary-time analysis as compared to 
the one based on centroid dynamics. 

To study the response to an externally applied field, 
a homogeneous force F is applied to each particle by 
adding a term — F^2 n X n to H in Eq. In Fig. [21 the 
FK chain's sliding velocity is shown as a function of the 
driving force for m — 0.02. This is the value of m, where 
the gap closes. It is found that the response in v is lin- 
ear with F at very small and very large F with different 
values for the effective damping coefficient 7 c ff — F/v. 
While (quantum) continuum approximations predict 7 e ff 
to be zero |7j, the chain's discreteness is known to change 
this property in classical systems, because kink-phonon 
interactions damp solitons |27j ■ Interestingly, 7 e g is inde- 
pendent of temperature T as T approaches zero in both 
linear regimes. This rules out the possibility that thermal 
fluctuations assist the system to overcome energy barriers 
at small F. Our observation is yet in partial contrast to 
the behavior of Luttinger liquids, which show v/F cx T a 
with a > 0. 

Another similarity with the dynamics of the classi- 
cal (discrete) FK model is seen at an external force 
F = 0.02, where the friction-velocity relation exhibits a 
cusp, separating a low-friction, low-velocity from a high- 
friction, high- velocity regime |2t|. Finally, at external 
forces F > 0.025, the response crosses over into a com- 
pletely unpinned regime, in which the motion of the FK 
chain is mainly opposed by the drag coefficient associated 
with the heat bath. In that regime, soliton-related dis- 
sipation mechanisms become unimportant. Transport of 
chains consisting of beads with masses m > m c requires 
either finite temperatures and/or finite forces. For ex- 
ample for m — 5to c , our numerical data provides an up- 
per bound for the mobility or inverse effective damping 
7-fjV = 0.1) < 2 • 10- 6 7 -ff (m = 0.02). 

The previous analysis shows that PIMD can distin- 
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FIG. 2: Sliding velocity v as a function of externally applied 
force F for mass m = 0.02, system sizes JV = 32 and N = 64, 
and temperatures T = 0.064 Vo/ks and T = 0.032, starting 
with a system at rest and slowly increasing F. The inset is 
an amplification of the small- u regime of the main figure. For 
F < 0.01, increasing or decreasing the thermostat (i.e., the 
externally imposed damping) by a factor of two has almost 
no effect on the v(F) curve. 



guish well between zero gap and finite gap. It is now 
possible to generalize the treatment for which the quan- 
tum mechanical ground state (of the continuum model) 
is not known analytically. One of the topics of current 
scientific interest is the interplay of disorder and quan- 
tum fluctuations 0|. Here we want to investigate the 
effect of disorder on the spectral properties of the quan- 
tum FK chain. This is done by replacing the external 
potential Vq cos(x n /b) in Eq. with a random poten- 
tial Uh{x), which is constructed in the following way: We 
add patches of the functional form Vq cos(x/b) of length 
nb where the underlying domain is chosen randomly to 
be either [0, nb] or [nb, 2nb). The patches are shifted by a 
constant so that no discontinuity in the potential occurs. 
This leads to {{V(x + Ax) - V(x)} 2 ) cx Ax 2H on scales 
b « Ax « nbN . The corresponding surface rough- 
ness exponent H takes the value H — 1/2. In order to 
construct a random potential with H = 0, the potential 
is either zero on a length of nb or - with same probabil- 
ity - it takes the functional form of Vb{l + cos(a;/6)} on 
the interval — n < x/b < n. Our numerical analysis of 
the results for the case H = suggests a finite gap for 



(H=0) 



and zero gap for m < m. 



(H=o) 



a mass m > m c „„ , . ... 

agreement with the predictions by Gorokhov et al. [15j. 
For the model potential under consideration, one phonon 
branch is observed at very large m, which then forms a 
broad band for masses larger than but in the order of m c . 
At m < m c , there is only one relatively narrow branch, 
which can be well described with Eq. and fco = 0. 

The situation is more complex for disordered poten- 
tials with roughness H = 1/2. In particular, we find that 
classical and quantum system behave qualitatively dif- 
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FIG. 3: Left: Zero wavevector spectrum for the quantum 
(m = 0.02, thick line) and classical (fcsT = 1.635 Vo, thin 
dashed line) FK model embedded in a disordered potential 
of roughness H = 1/2. T was chosen such that the average 
kinetic energy was the same in the classical as in the quantum 
case. Right: Full spectrum of the velocity correlation function 
for the quantum system. Average over 12 different disorder 
realizations. 



ferent. In agreement with the predictions by Chauve et 
al. there is no indication for the classical mobility 
to become zero at finite temperatures. In contrast, the 
quantum mechanical mobility always appears to be zero 
for H = 1/2, no matter how small m is. Whenever, we 
were prompted to believe the quantum H = 1/2 system 
was depinned, doubling the system size dramatically re- 
duced the velocity- velocity spectrum C(q, u>) near zero 
wavelengths and zero frequencies. The zero- wavevector 
spectrum C(q — 0, lu) for classical and quantum mechan- 
ical FK chains is shown in the left-hand side of Fig. [3] 
where the qualitative difference of C{q, to) in the limit 
(q,td) — > (0,0) is bourne out. (Note that due to statis- 
tical uncertainties and finite system size, one can never 
observe C(0,0) = in a computer simulation.) 

The spectrum at finite wavelengths of the quantum 
system, which is shown on the right-hand-side of Fig. [3] 
shows two phonon branches. The figure is the average 
over 12 different disorder realizations. A single realiza- 
tion shows much sharper lines. It appears that the system 
is not self-averageing. 

In conclusion, we find that the regular zero- 
temperature, discrete, quantum sine Gordon model has 
an insulating region with m > m c and finite gap, while 
for m < m c the gap closes resulting in finite resistance. 
Similar behaviour is found for a quantum chain on a dis- 
ordered potential with roughness exponent H = 0. How- 
ever, H = 1/2 automatically leads to a zero-temperature 
insulator, even though small thermal fluctuations are suf- 
ficient to induce creep. 
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